Scaling for the Percolation Backbone 
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We study the backbone connecting two given sites of a two-dimensional lattice separated by an 
arbitrary distance r in a system of size L. We find a scaling form for the average backbone mass: 
(Mb) ~ L d ° G(r / L) , where G can be well approximated by a power law for < x < 1: G(x) ~ x^ 



with V = 0.37 ± 0.02. This result implies that (M B ) ~ £ 



for the entire range < r < L. 



We also propose a scaling form for the probability distribution P(Mb) of backbone mass for a given 
r. For r &s L, P(Mb) is peaked around L dB , whereas for r <^ L, P(M B ) decreases as a power law, 
Mg TB , with Tg ~ 1.20 ± 0.03. The exponents tp and tb satisfy the relation ip — d B (r B — 1), and ip 
is the codimension of the backbone, ip — d — d B . 
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PACS numbers: 64.60.Ak, 05.45. Df 



The percolation problem is a classical model of phase 
transitions, as well as a useful model for describing con- 
nectivity phenomena, and in particular for describing 
porous media At the percolation threshold p c , the 

mass of the largest cluster scales with the system size L 
as M ~ L d i . The fractal dimension df is related to the 
space dimension d and to the order parameter and cor- 
relation length exponents /3 and v by df = d — fl/v 0-[jJ. 
In two dimensions, df = 91/48 is known exactly. 

An interesting subset of the percolation cluster is the 
backbone which is obtained by removing the non-current 
carrying bonds from the percolation cluster 0]. The 
structure of the backbone consists of blobs and links 
0.0 0. The backbone can in fact be further partitioned 
into subsets according to the magnitude of the electric 
current carried HPL The backbone is relevant to trans- 
port properties |lf||] and fracture . The fractal di- 
mension ds of the backbone can be defined via its typ- 
ical mass Mb, which scales with the system size L as 
M B ~ L dB . The backbone dimension is an independent 
exponent and its exact value is not known. A current 
numerical estimate Jl|] is d B = 1.6432 ± 0.0008. 

The operational definition of the backbone has an in- 
teresting history Jl|-§]. Customarily, one defines the 
backbone using parallel bars, and looks for the percola- 
tion cluster (and the backbone) which connects the two 
sides of the syste m H j. A different situation arises in oil 
field applications |]l2| , where one studies the backbone 
connecting two wells separated by an arbitrary distance 
r. This situation is important for transport properties, 
since in oil recovery one injects water at one point and 
recovers oil at another point |]l2|| . From a fundamen- 
tal point of view, it is important to understand how the 
percolation properties depend on different boundary con- 



ditions. 

We study here the backbone connecting two points sep- 
arated by an arbitrary distance r in a two-dimensional 
system of linear size L. One goal is to understand the 
distribution of the backbone mass Mb (r, L) , and how its 
average value scales with r and L in the entire range 
< r < L. 

We choose two sites A and B belonging to the infinite 
percolating cluster on a two-dimensional square lattice 
(the fraction of bonds is p = p c = 1/2). A and B are sep- 
arated by a distance r and symetrically located between 
the boundaries Using the burning algorithm, we 

determine the backbone connecting these two points for 
values of L ranging from 100 to 1000. For each value of L, 
we consider a sequence of values of r with 2 < r < L — 2. 
In order to test the universality of the exponents, we per- 
form our study on three lattices: square, honeycomb and 
triangular lattice. For simplicity, we restrict our discus- 
sion here to the square lattice, as we find similar results 
for the other two lattices. 

We begin by studying the backbone mass probability 
distribution P(Mb)- We show that P(Mg) obeys a sim- 
ple scaling form in the entire range of r/L, 



(1) 



where F(x) is a scaling function, whose shape depends 
on the ratio r/L. 

For r ~ L, it seems reasonable to assume that P(Mb) 
will be peaked around its average value < Mb >^ L dB . 
The data collapse predicted by Eq. (0) is represented in 
Fig. 1(a). In this case, the scaling function F is peaked 
at approximately L dB . 
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However, the case r <C L is far less clear. In fact, we ex- 
pect for r <C L that the backbone mass fluctuates greatly 
from one realization to another, since its minimum value 
can be r and its maximum can be of order L df . Fig. 1(b) 
shows a log- log plot of P(Mb). It has a lower cut-off of 
order r (since the backbone must connect points A and 
B) and a upper cut-off of order L dB . We find good data 
collapse (Fig. 1(c)), which indicates that the scaling func- 
tion F is a power law in the range from r dB to L dB , with 
exponent approximately tb — 1.20 ±0.03 (there is a cut- 
off at Mb ~ L ds not shown here). The exponent tb is 
connected to the blob size distribution ||] since typically, 
the two sites belong to the same blob, and the sampling 
of backbones is equivalent to sampling of the blobs. From 



which means that tp is the codimension of the fractal 
backbone. 

To summarize, we find that for any value of r/L, the 
scaling form, Eq. (Q), for the probability distribution is 
valid. The shape of the scaling function F depends on 
r/L, being a peaked distribution for r ss L, and a power 
law for r <C L. The average backbone mass varies with 
r and L according to Eq. (Eh. For fixed system size, it 



varies as (Mb) 



(for < r < L) . The value of ip is 



small (ip f» 0.37) indicating that the backbone mass does 
not change drastically as r changes. On the other hand, 
the exponent governing the variation of (Mb) with L for 
fixed r is expected to be larger, with (Mb) ~ r ds ~^ . 
This exponent cLb — ip is not equal to the fractal dimen- 
sion cLb of the backbone, but is smaller by an amount 
equal to ip. 



TB- 



(2) 



This relation gives the estimate tb — 1-22 in good agree- 
ment with our numerical simulation. 

We note that for larger values of Mb, a "bump" (indi- 
cated by an arrow on Fig. 1(b)) located at approximately 
L ds appears and assumes increasing importance when r 
approaches L. 

We now study the average backbone mass (Mb)- From 
dimensional considerations, the r dependence can only be 
a function of r/L. We thus propose the following Ansatz: 



(M B (r,L))=L dB G(^) 



(3) 



In Fig. 2(a), we show a double logarithmic scale Mb ver- 
sus r for different values of L. In order to test the Eq. (||), 
we scale the data of Fig. 2(a). The data collapse is ob- 
tained using cLb — 1-65 and is shown on Fig. 2(b). This 
(log-log) plot supports the scaling Ansatz (||). Moreover, 
one can see that the scaling function G is, surprisingly, a 
pure power law on the entire range [0, 1], with exponent 
ip = 0.37 ±0.02. 

The results (Q) and (|J) are consistent, since if (|l|) 
holds with a power law behavior for the scaling func- 
tion F(x) ~ x~ TB for x > 1, and F(x) = for x < 1, 
then the average mass is given by 



(M B (r,L)) 



L B 



F 



M\ dM 



-M. 



(4) 



Assuming that L/r is large enough, the integral in (ji|) 
can be approximated as L dB -^r^, where 



ip = cLb{tb - 1) 



(5) 



In our simulation tb ~ 1-20 ± 0.03, which leads to the 
value ip w 0.33 ± 0.05 in reasonable agreement with the 
value measured directly on the average mass. 

Moreover, using Eq. @ together with Eq. (||), we ob- 
tain 
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FIG. 1. (a) Data collapse of P(M B ) using Eq. (1) for three 
different value of r ~ L. (b) Probability distribution of the 
backbone mass for L = 1000 and r = 2 (computed with 10 s 
configurations). The exponent tb is obtained by a linear fit 
over the range 30 < Mb < 3 x 10 4 and the error bar on tb 
is around 0.03. The arrow denotes the fact that Mb peaks as 
L dB . (c) Data collapse of P(M B ) for L = 400 using Eq. (1) 
for three different values of r. 
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FIG. 2. (a) Log-log plot of the average backbone mass 
{M B ) versus r for four different values of L. (b) Data from 
Fig. 2(a) collapsed with the use of the scaling form proposed 
in Eq. (2). The error on ip is typically 0.02. 
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